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Approximation only by derivative (or more generally momentum) expansions, com- 
bined with reparametrization invariance, turns the continuous renormalization 
group into a set of partial differential equations which at fixed points become 
non-linear eigenvalue equations for the anomalous scaling dimension r\. We review 
how these equations provide a powerful and robust means of discovering and ap- 
proximating non-perturbative continuum limits. Gauge fields are briefly discussed. 
Particular emphasis is placed on the role of reparametrization invariance, and the 
convergence of the derivative expansion is addressed. 



This talk is about derivative (or more general momentum) expansions 
as applied to the renormalization group, in a quantum field theory setting. 
The motivation is simply this: I want to construct analytic approximation 
methods with as much reliability and accuracy as possible even when there 
are no obviously small parameters, e.g. e = 4 — D, l/N etc. with which one 
could expand perturbatively. (Here, D is space-time dimension, and N is 
the number of components of the field.) In other words, I want to look for 
approximations that work in a genuinely non-perturbative setting. Now, as 
Wilson was instrumental in demonstrating, the continuum limit of a quantum 
field theory can non-perturbatively be best understood in terms of the flow of 
aii. effective action 5a [<p] under lowering an effective U.V. (ultra-violet) cutoff 
Alii Thus in this framework one works with a flow equation that generically 
takes the form 

A^aM=^[5a] , (1) 

where the cutoff is implemented through some effective U.V. cutoff function 
Cj/y(g, A). (Here q stands for momentum, and the above equation is referred 
to as the continuous, or momentum space, renormalization group.) Scale in- 
variant continuum limits (thus massless field theories) are then simply given 
by fixed points: [j] 

A^5aM = . (2) 

The massive continuum limits follow from tuning the relevant perturbations 
around these fixed points. In such a setting one realises that various approxi- 

a once all quantities have been rewritten in terms of dimensionless quantities, using A 
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mations can be made quite easily that preserve the structure of the continuum 
limit, while in other frameworks (for example when using truncations of Dyson- 
Schwinger equations) the continuum limit, equivalently renormalisability, is 
almost inevitably destroyed^ 

What are the possible approximations? The first thought is to try trun- 
cating the space of interactions to just a few operators, however this results 
in a truncated expansion in powers of the field tp (about some point). Such 
an approximation can only be sensible if the field ip does not fluctuate very 
much, which is the same as saying that it is close to mean field, i.e. in a setting 
in which weak coupling perturbation theory is anyway valid. Studying the 
behaviour of truncations in a truly non-perturbative situation, one finds that 
higher orders cease to converge and thus yield limited accuracy, while there is 
also no reliability - even qualitatively - since many spurious fixed points are 
generated^ 

This situation should be contrasted with truncations to a few operators 
in the real space renormalization group of spin systems, such as block spin 
renormalization group of the Ising model. Such-truncations were extensively 
studied in the pastp and could be very accurately. A modern variant produces 
spectacularly accurate results in. low dimensionsQ The powerful Monte-Carlo 
renormalization group methods^ are also based on such truncations. In the 
case of such simple discrete systems however, the expansion effectively results 
in a short distance expansion of the effective action, for example by keeping 
only the finite number of interactions linking nearest neighbours, then next- 
to-nearest neighbours, and so on. 

The analogous expansion in our continuum case is, for smooth cutoff 
functions Cjjv, a derivative expansion of-J^A [v 3 ] (equivalent to a Taylor ex- 



CuvisL, A) = 9(A — q) it is an expansion in momentum scale ~ |<9/A| = p/A, 
where the coefficients are not analytic in p^ but rather, non-trivial func- 
tions of the angles between various momenta which must be determined self- 
consistently through the flow equationEJ [The non-analyticity is a purely tech- 
nical problem that is induced by the non-analyticity of 8(A ~ q).] At any rate, 
in both cases such a short distance expansion - where no other approximation 
is made - seems a particularly natural approximation to try, and in view of the 
discussion above, sensible results might well be expected providing the Wilson 
effective action S\ is 'sufficiently well behaved': thus the approximation would 
fail if the higher derivative terms are not in some sense small, but this would 
indicate that a description in terms of the given field content is probably itself 
inappropriate and other degrees of freedom should be introduced. This is an 
important point, to which I will return later. 



pansion in the momenta 




cutoff functions 
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Consider the case of O(N) invariant scalar field theory, the so-called N- 
vector model. I shall start by using a sharp cutoff and making the simplest 
such approximation - keeping only a potential interaction: 

S A ~Jd D x {\{d^ a f + V{p,K)} . (3) 

After appropriate scalings to dimensionless combinations, the flow equation is 
found to be: 

d 

- A—V + d(pV - DV = ln(l + V ) + (N - 1) lnfl + V lw) . (4) 

OA 

Here (p = y/ip a tp a , I = d/dtp, and I have introduced the dimension of the field 
ip: d = D/2 — 1. (Since we have thrown away all momentum dependent cor- 
rections, r) = in this approximation.) The N — oo case of this equation was 
already derived by Wagner and Houghton in their paper introducing the sharp- 
cutoff flow equation!^ and subsequently the general N case was..proposed as 
a "local potential" approximation by_Nicoll, Chang and StanleyEj It has since 
been rediscovered by many authors t3 especially Hasenfratz and HasenfratzEij 
Nevertheless this equation, and its smooth cutoff sisters, have a number of 
beautiful properties that have not been pointed out by previous workers. 
First note that the fixed point equation for V((p, A) = V(<p), 

dtpV - DV = ln(f + V") + {N- l)ln(l + V'/ip) , (5) 

determines by itself at most a countable set [] of sensible fixed point potentials, 
each of which can be identified with approximations to the exact fixed points. 
This is not obvious because eqn.(||) is a second order ODE (Ordinary Differen- 
tial Equation) and therefore in some neighbourhood of some starting value tp 
we can construct a continuously infinite two parameter set of solutions. Actu- 
ally, all but a countable number of those solutions are singular! For example, 
consider the case D — 4 and N — 4 (the Higgs field in the Standard Model). 
Obviously from (§), V'(0) = is necessary if the potential is not to be singular 
at ip = 0. We can choose a value of V(0) and then (numerically) integrate out 
to ip > using eqn.(^). One discovers almost always that at some point (p = ip c 
a singularity in V is encountered, after which the potential ceases exist (or at 
least is complex for tp > tp c ). The first graph in fig.l is a plot of p c against 
V(0). We see that only the (trivial) Gaussian fixed point solution V(tp) = 
exists for all values of the field. If the same is done for the case D = 3 and 

b D = 2 is an exception: see later 
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N — 1, we get the second graph in fig.l. In this case there is also one non- 
trivial non-singular solution, corresponding to the famous Wilson-Fisher fixed 
point (Ising model universality class). I have checked all cases D = 3,4 and 
N = 1, 2, 3, 4; they reproduce the standard fixed pointsH 




Figure 1: Plots of tp c against V(0) for D = N = 4, and D = 3 and JV = 1. 



Studying eqn.(|g), one can convince oneself that the only way that V can 
satisfy this equation as (f — ► oo is if V(<p) ~ ip D l d . Together with V(0) = 0, we 
now have two boundary conditions and thus we should expect only a countable 
number of solutions from the second order ODE. For D = 2, exceptions arise 
(due to d = 0), thus for N = 1 one obtains a semi- infinite continuous line of 
fixed points with periodic potentials, but these may be identified with critical 
sine-Gordon modelsEI On the other hand, there is no need to impose V^'(0) = 
when N = 1. This just allows a constant phase shift on critical sine-Gordon 
potentials, while in higher than two dimensions the power law constraint on 
V now holds separately for tp — > oo and tp — > — oo (i.e. with possibly different 
coefficients). Nevertheless, we have confirmed that in this larger space, there 
is still only the one non-trivial fixed point in three dimensions. 

While all this just reproduces the standard lore, note nevertheless how 
powerful the method is: the infinite dimensional space of all possible potentials 
V(ip) has been searched for continuum limits. Clearly this is much more than is 
possible with other methods! Also, the continuum is actually accessed directly 
without the need to go through the construction of introducing an overall cutoff 
Ao, a bare action S\ , and then taking the continuum limit Ao — * oo. 

These properties are true also when the approximation is applied to the 
massive theory. In this case one must determine the form of the perturbations 
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about the fixed point. One can write 

V(<p,A)=V(<p)+ev x (<p)A- x , (6) 

where eqn.(^j) is linearised in e and separation of variables has been used. Now, 
v\(tp) satisfies a linear second order ODE, and once again we appear to have a 
continuously infinite set of solutions, and for all choices of A. In this case the 
crucial observation is that if, beyond the linearised level, the scale dependence 
of the perturbation is to be absorbed into an associated coupling g\ (A), i.e. 
a renormalised coupling corresponding to universal self-similar flow abopi the 
fixed point, then v\(ip) has to behave as V\(ip) ~ (p( D ~ x )/ d as ip — > oo. EJ For 
the same reasons as before, this typically allows only a countable number of 
solutions, but this time we also have linearity, which implies a normalization 
condition can be set, overconstraining the equations and resulting in quanti- 
zation of A. (Again, d — provides an exception - resulting in more general 
perturbations with e.g. exponential or periodic behaviour. It is worth remark- 
ing that for d ^ 0, the power laws given above are the unique powers required 
so that the physical [] v\ and V are independent of A, and therefore obtain a 

Son-trivial finite limit as A — > 0; this limit gives the Legendre effective potential 
and thus the equation of state.) 

Consider now the derivative expansion at 0(d 2 ). In this case we need to 
use a smooth cutoff (as already discussed) . The effective action takes the form 

S A ~ JdPx + \{d^ a fK{ v ) + \{ v a d^ a fZ{ v )} (7) 

(where this last term is required only for N ^ 1). In this case the fixed point 
equations are a set of coupled second-order non-linear ODEs, one for each coef- 
ficient function (V, K and Z). This pattern holds to all orders of the derivative 
expansion. As previously, one can argue for specific power law behaviours for 
V, K and Z , and that typically only a countable number of non-singular so- 
lutions exist. But now there is another parameter to determine: the critical 
exponent 77 from the anomalous scalina-of-^. The exact renormalization group 
has a reparametrization invariance reflecting the fact that physics is 
independent of the normalization of the field, and this extra invariance turns 
the fixed point equations into non-linear eigenvalue equations for 77, because it 
allows a normalization condition (e.g. K(0) = 1) and thus quantization of r\ (in 
a similar way to the linear case above for perturbations.) There is a problem 
however: the derivative expansion generally breaks the reparametrization in- 
variance, with the result that 77, v, to etc. depend on some unphysical parameter 
such as K(0), the normalization of the kinetic term. Ell 

c i.e. in the original dimcnsionful variables. 
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Consider the Polchinski forma of the Wilson flow equation. Schematically, 
dS A ldA uv (6S A 6S A S 2 S A , ! ,5S A \ 

where Auv(q,A) = Cuv/q 2 , and the total action S a, ~ | ip.A^y.ip +S\ nt . 



This equation is simaly related to the Wilson equation^ through tp i— > VCuv ¥> 



and H = —S\ , EHIj but in contrast to Wilsons, it has the intuitively nice 
property that eigen perturbations about the Gaussian fixed point are precisely 
polynomials in the field and its derivatives. Now in general, the repararo£jtriza- 
tion symmetry is given by a complicated functional integral transformpl so it 
is not surprising that a truncated derivative expansion destroys it - and in fact 
it is far from clear how to approximate at all in a way which preserves it. 

(Let me emphasise that with broken reparametrization invariance, defining 
z(ip) through SjC* ~ / V(<p) + \z((p){d^(p a ) 2 +■■■, the results depend on the 
value z(0). This dependence on z(0) has sot been recognized by authors who 
set z(0) = with insufficient justificationE3) 

Fortunately, for two special forms of cutoff function, reparametrization 
invariance may be linearly realized. These correspond to either sharp cutoff 
Cuv = 0(A - q) or a power law smooth cutoff Cuv ~ 1/[1 + (q/A) 2K+2 ], 
where k is some non-negative integer. (The underlying reason is that such 
cutoffs are left invariant by a subgroup of linearly ptealised 'universality sym- 
metries' that map between different schemes in (||)P3 but it would take us too 
far afield to explain this.) However, a direct derivative/momentum-scale ex- 
pansion of the_Wjlson flow equation leads to singular coefficients with both of 
these cutoffs crcl't3'E3 The way out of this difficulty is to recognize that, from 
the form of the. right hand side of eqn. (||) , the Wilson effective action has a 
tree structures It is the.JTay.lnr expansion of the corresponding propagators 
that causes the problcm&BE3c!J Therefore, to overcome this difficulty we first 
pull out the one particle irreducible parts. 

It can be showncl that the one particle irreducible parts of S A are generated 
by a Legendre effective action T A [<p] equipped with infrared cutoff Cm = 
1 — Cjjvj *- e - related in the usual way to a partition function except that the 
bare action is modified to S Ao = ^tp.Ajl.(p where Air = Cm/q 2 . 



Whence, the 'Legendre flow equation' follows; 
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(9) 



In this fopii, the sharp cutoff limit enjoys the simplest reparametrization 
invariancet3 <p i— > aip. For the smooth power law cutoff, certain momeiu 
turn independent linear transformations on other quantities are also requiredB 
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Since these transformations are linear and momentum independent, a trun- 
cated derivative expansion now preserves the reparametrization invariance, 
while depyative/momentum-scale expansion of T™* is well-defined with these 
cutoffspE3 because it results in Taylor expansion of the self-energy, rather than 
the propagator itself. 

Returning now to 0(<9 2 ), we need the smooth cutoff and choose the integer. 
k as small as possible, to maximise the accuracy of the derivative expansion^ 
Table 1 displays the results obtained in D = 3 dimensions@'Ej for r\ and ^,[] 
and for the first correction to scaling exponent u>. Only the expected fixed 
points are found. For N = 1 • • • 4, the results are fair for the local potential 
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Table 1: Critical exponents of the three-dimensional Wilson-Fisher fixed point. The first 
two ordc«3-|(j if i th e derivative expansion are compared to a combination of the worlds best 
cstimatcsM'EjE] with their errors, where available. r\ is identically zero for all N at O(d ). 

approximation and already quite good at 0(d 2 ), with the exception of u> which 
gets worse at 0(d 2 ) when N > 1. For the large N cases however, ui is better 
estimated, v is not improved at 0(d 2 ), while n is dramatically underestimated 
- eventually by about a factor of 10. Something is going wrong particularly at 
large N. Actually, the results for N — oo are guaranteed to be right because 
in this case the O(d ) approximation is exacttJ'E2l so the problem only appears 
in the approach to this limit. I believe that this is an example of a case 
where all the appropriate fields have not been included in the effective action. 
Indeed, it is known that at large N, a massless bound state field also exists 

"from which all other exponents (excepting correction to scaling exponents) follow, because 
all hyper-scaling relations are here exactly preserved 
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at the critical point oE3 We should thus expect that a derivative expansion 
is ill-behaved, because the vertex functions are hiding within them the effects 
of this integrated out massless field. To ameliorate this behaviour, we should 
include the bound state explicitly as an O(N) singlet field, then amongst the 
new set of fixed points in this enlarged space will be one with the same universal 
properties as the original N vector model, but with better behaved derivative 
expansion properties. Similar considerations should apply to fixed points with 
fcrmions, particularly since the bound state fields here also correspond to the 
order parameter (a.k.a. fermion condensate) EZI 

The most impressive example so far however, is provided biL-the case of 
D = 2 dimensions. It had been conjectured by Zamalodchikov Ej that there 
should exist an infinite series of multicritical points in the two dimensional 
Z2 (i.e. if <-> — cp) symmetric theory of a single scalar field, corresponding to 
the so-called unitary minimal models in CFT (conformal field theory). How- 
ever a verification of this conjecture is in practice well outside the capabilities 
of the standard non-perturbative approximation methods: the corresponding 
e expansions are so badly behaved as to be useless,^] with similar difficulties 
expected in resummed weak coupling perturbation theory, while lattice meth- 
ods suffer from difficulties locating and accurately computing the multicritical 
points in these high dimensional bare coupling constant spaces.^] All these 
methods also get rapidly worse with increasing operator dimension. In con- 
strast, at 0(d 2 ), the lowest order at which a fixed point search through 77 > 
can be done, we uncover the multicritical points, and only these^J and find 
an agreement with CFT that improves with increasing multicriticality and 
dimensionP The results are displayed in table 2, for the first 10 (multi)critical 
points and up to the first 10 operators. We see that there is a remarkable 
agreement between these thus lowest order results and CFT, spanning over 
two orders of magnitude. The worst determined number is rj for the tricrit- 
ical point, which is only accurate to 33%, but this gradually improves as to 
increases. (At m = 11, rj is off by 21%.) v is worst determined at to = 3 
(13%) after which all are determined to error less than 2% and decreasing 
with increasing to. Indeed, the best determined number is v for m = 11, 
which is accurate to 0.2%. The worst determined operator dimension is the 
3 rd at to = 5 (25%), after which errors decrease with increasing to and/or 
increasing dimension, so that all the rest have errors in the range 9% - 22%. 



e E.g. already with the tricritical point, 0(e 2 ) underestimatep-p by a factor ~ 1/100. The 
situation gets factorially worse as multicriticality is increasedE3 

Indeed to date, only the lattice computation of the two lowest operator dimensions around 
the tricritical point, has been attemntedrl 

3 The search was restricted to real ^—symmetric V and K, with K(Q) > 0. 
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Table 2: Critical exponents, and dimensions and Z2 parities of the 10 lowest dimension 
operators, for the first 10 multicritical points, m = 2 ■ ■ ■ 11. Thus, m = 2 is Ising, m = 3 is 
tricritical Ising, m = 4 is tetracritical (a.k.a. three-state Potts model), etc. For each m, the 
0(d 2 ) answer is in the first row and the associated exact CFT result is in the second row. 



Note that the low dimension operators at high order of multicriticality corre- 
spond to renormalization group eigenvalues which agree with CFT to better 
than 3 significant figures. Fig. 2 shows the fixed point solutions for the first 
three critical points. 

We also found some irrelevant operators which cannot be matched to op- 
erators in the CFT minimal models. This is the reason for the blank spaces in 
some corresponding CFT parts of the table. Examples of such operators have 
also been found (in the correction to scaling at the Ising critical point) by e 
expansion and fixed dimension resummed perturbation theorypj and argued 
for in exact treatmentsEIJ Further work is required to understand their true 
significance. 

Finally, what are the problems with this approach? The treatment of gauge 
theory present special problems: the formulation of a flow equation which 
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Figure 2: The potential V((f), and kinetic term coefficient function K(ip), for the two di- 
mensional Ising, tricritical Ising, and tetracritical Ising model fixed points. 



properly treats the quantum aspects of gauge invariance non-perturbatively, 
and secondly, the construction of reliable approximations. The solution to the 
former must proceed in one of two ways: by allowing the cutoff to break the 



gauge invariance and then a 
(via broken Ward identities) 
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mpting to recover it as the cutoff is removed 
_ or by generalizing the flow equations so that 
gauge invariance is not broken by the cutofftS Whilc-the first method can be 
shown to work— to all orders in perturbation theorypa it seems hopeless non- 
perturbativelya ft is surely through an exact preservation of the quantum 
gauge invariance_that real progress will be made; a useable generalization of the 
pure U(i) caseE3 may well be possible. Note however, that for the interesting 
cases (non-Abelian and greater than two dimensions), derivative expansions 
per se are anyway impractical due to the large number of independent gauge 
invariant combinations that can be formed even at lowest non-trivial ordered 
Instead, much more appropriate forms of approximation deserve study in this 
case, such as large N methods. 

Returning to non-gauge theories, apart from the practical problem that 
higher orders in the derivative expansion get rapidly more complicated, one 
problem that this method shares with all other approximations to the renor- 
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malization group is an unphysical dependence on the choice of cutoff func- 
tion. This dependence is not really a problem however if instead it is used 
to estimate a rough lower bound on the error of the approximation, and thus 
test the numerical reliabilityE] For example, we can compare the O(d ) N = 1 
results ixLtable 1, to the corresponding results for sharp cutoff, v — .70 and 
lo = .60E3B Another problem that has already been mentioned is that gener- 
ically derivative expansions depend also on one unphysical parameter in the 
fixed point solutions, due to loss of reparametrization invariance. In a similar 
way this is not a problem (if one knows about it 1 ) but rather can be used to test 
numerical reliability in such approximationslaO The real question that needs 
to be answered ultimately, is whether the derivative expansion exists to all or- 
ders and converges, because of course if it does, these problems have to lessen 
and eventually disappear as the expansion is pushed to higher orders. This 
question is hard to answer in generality, but rather straightforward to analyse 
perturbatively for a specific theory. In table 3, I show which results converge 
for the (3 function of the O(N) invariant scalar field theory in D = 4 dimexbj 
sions, computed to two loops with different forms of derivative expansiont3 
Also summarised are which expansions preserve reparametrization invariance. 



Variant. 


Rcpar. Inv. 


Convergence. 






one loop 


two loops 


Wilson/Polchinski 
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X 


X 


Legendre Power 
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X 
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X 


V 


V 


Legendre Sharp 


V 


V 


V 



Table 3: Properties of derivative expansions for different forms of flow equation and cutoff: 
Wilson or Polchinski flow equation with any smooth cutoff, Legendre flow equation with 
the Power law cutoff described in the text, Legendre flow equation with some Faster falling 
cutoffs, and momentum scale expansion of the Legendre flow equation with Sharp cutoff. 

A direct derivative expansion of the Wilson effective action does not con- 
verge already at one loop: it is again a result of expansion of the propagators 
inside the effective action (mentioned earlier). On the other hand, the Legen- 
dre flow equations give the exact answer to the one loop j3 function already 
at O(d )a Only the last two methods however converge at twq-Jpops. In par- 
ticular, the sharp cutoff case converges very fast at two loops^ and since it 

h and analogously perturbation theory to a certain extent, through scheme dependence. 
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also preserves a simple reparametrization invarianr:e, further work to overcome 
the practical difficulties in its implementation E3 certainly seems called for. 
Although convergence even to all orders in perturbation theory, is no guar- 
antee of convergence non-perturbatively (truncations in increasing powers of 
the field trivially must converge at any fixed order of perturbation theory, but 
as we have seen, fail to do so non-perturbatively), the natural conjecture is 
that these last two methods do converge non-perturbatively. The fact that the 
Legendre O(d ) approximation is exact at N = oolldo lends further support 
to this conjecture.| 

Of course, negative answers to convergence do not exclude the first two 
methods from being good model approximations at low orders. Although I 
have concentrated on the Legendre flow equation with power law cutoff, deriva- 
tive expansions of the Wilson / Polchinski equation are distinguished by their 
relative simplicity. Perhaps, by generalising the derivative expansion, one can 
preserve this simplicity while also preserving more of the structure of the exact 
renormalization group. 

We remind that a large number of references to otheij— work on the low- 
est order sharp cutoff approximation have been collected I 1 Pl.We collect here 
corresponding smooth cutoff versions not so far mentionedpE3 and similarly 
attempts to go beyond leading order in the derivative expansionEZI There are 
also a number of examples that entertain the idea of derivative expansion but 
in practice make further truncations. 

In conclusion, the derivative (momentum scale) expansion methods - where 
no other approximation is made - are potentially very powerful, particularly in 
genuinely non-perturbative settings where all other methods fail. In contrast 
to more severe truncations, all these variants are robust, in the sense that no 
spurious solutions have been found, while especially the Legendre flow equa- 
tion with power law cutoff yields very satisfactory numerical accuracy at low 
orders. The full potential of these methods is by no means yet realised, and 
much more theoretical progress on their properties is possible and expected. 
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l E.g. the same is not true of truncations in powers of the field, to any order. 
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